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Abstract 



It is pointed out that the enhancement of the solar neutrino rate in a real 
time detector like Super-Kamioka, SNO or Borexino due to neutrino oscilla- 
tions in the moon during a partial or total solar eclipse may be observable. 
The enhancement is calculated as a function of the neutrino parameters in 
the case of three flavor mixing. This enhancement if seen, can further help to 
determine the neutrino parameters. 
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I. INTRODUCTION 



The sun is a copious source of neutrinos with a wide spectrum of energies and these neu- 
trinos have been detected by terrestrial neutrino detectors, although at a rate lower than ex- 
pected from theoretical calculations. A new generation of detectors with high counting 
rates will soon be producing abundant data on solar neutrinos. Mixing and the consequent 
oscillations among the neutrinos of different flavors is generally believed to be the cause of 
the reduced intensity of neutrino flux detected on the earth. However, neutrino-oscillation 
is a complex phenomenon depending on many unknown parameters (six parameters for 
three flavors u e , v^, u T ) and considerable amount of experimental work and ingenuity will 
be required before the neutrino problem is solved. 

Hence it would be desirable, if apart from direct detection, we can subject the solar 
neutrino beam to further tests by passing it through different amounts of matter, in our 
attempts to learn more about the neutrinos. Nature has fortunately provided us with such 
opportunities: (1) Neutrinos detected at night pass through the earth. (2) Neutrinos de- 
tected during a solar eclipse pass through the moon. (3) Neutrinos detected at the far side 
of the earth during a solar eclipse pass through the moon and the earth. We shall call this 
scenario (3) a double eclipse. The scenario (1) has been studied in the literature rather 
extensively The purpose of the present work is to examine the scenarios (2), and (3). 

Two previous works ||[7| have discussed scenario (2), however both are incomplete in many 
respects. 

The plan of the paper is as follows. The relevant astronomy is presented in Sec.[H|. In 
Sec.[Hj we give the theory of the passage of the solar neutrinos through the moon and the 
earth, taking into account properly the non-adiabatic transitions occuring at the boundaries 
of the moon and the earth. In Sec.[IV| we present the numerical calculations of the neutrino 
detection rates during the single and double eclipses and the results. Sec.|V| is devoted to 
discussion. 
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II. ECLIPSES AND DOUBLE ECLIPSES 



Solar neutrinos are produced within the solar core whose radius is of order 1/10 of the 
solar radius and we shall approximate this by a point at the centre of the sun. What is 
required for our purpose is that the lunar disc must cover this point at the centre of the 
sun and so as far as the neutrino radiation is concerned , the solar eclipse is more like an 
occultation of a star or a planet by the moon. 

Astronomers characterize the solar eclipse by the optical coverage C which is defined 
as the ratio of the area of the solar disc covered by the lunar disc to the total area of the 
solar disc. For neutrino physics we require the distance du travelled by the neutrino inside 
the moon. Defining the fraction x = ^ where Rm is the lunar radius, x can be given in 
terms of C by the following formulae: 



eqn(|l|). The relationship between x and C so obtained is plotted in FigQ. When the lunar 
disc passes through the centre of the sun, C is 0.39 and the neutrino eclipse starts at this 
value of C.When the optical coverage increases above 39%, x rises sharply from zero and 
reaches 0.6 and 0.95 for optical coverage of 50% and 80% respectively. 

For any point of observation of the usual solar eclipse (which we shall call single eclipse) 
there is a corresponding point on the other side of the earth where a double eclipse occurs. 
With the coordinates labeled as (latitude, longitude), the single eclipse point (a, (3) is related 
to the double eclipse point (A, a) by the relations (see Fig.Q): 




(1) 




(2) 



Eqn(|2]) can be inverted to get the parameter z for a given C and this z can be substituted in 



A = a + 25 



a = 7i - 2 Q UT - (3 



(3) 



where, 
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2nf 

sin 5 = sin 23.5° sin( — ), (4) 
Ty 

Ty is the length of the year , zero of time t is chosen at midnight of autumnal equinox 
i.e. Sept. 21, and Qjjt is the angle corresponding to the Universal Time - UT. 

During a double eclipse, the neutrinos travel through the earth in addition. The distance 
(1e travelled by the neutrino inside the earth along the chord between the points (a, (3) and 
(A, a) as a function of time t, is given by 

2nt 

cIe = 2i?£:(sin A sin 5 + cos Acosc*) cos(— — )) (5) 

Td 

where Re is the radius of the earth and To is the length of the day. This is the same distance 
that needs to be calculated in the study of the day-night effect and a plot of this distance 
as a function of t is given in our earlier paper ||. 

Present and upcoming high statistics neutrino detectors expect to collect a few solar 
neutrino events every hour. As discussed in Sees. |111 A| and |ill B| , single and double eclipse 



can lead to enhancements of rates by upto two and a half times. Even with such large 
enhancements during the eclipse the signal may not exceed statistical errors, since each 
solar eclipse lasts only for a few hours. However they occur fairly often. As many as 32 
solar eclipses are listed to occur during the 14 year period 1996 through 2010. Global maps 
and charts are available || for location and duration of both the umbral and penumbral 
coverage. 

The eclipses during the 2 year period 1997 through 1999 are the following: 
1997 March 9 - Total Solar Eclipse 

1997 Oct 12 - Partial Solar Eclipse 

1998 Feb 26 - Total Solar Eclipse 

1998 Aug 22 - Annular Solar Eclipse 

1999 Feb 16 Annular Solar Eclipse 



1999 Aug 11 Total Solar Eclipse 

We have analyzed the five total/annular eclipses. We have not examined the partial solar 
eclipse of 1997 October 12, since we have not so far been able to procure the data for this 
eclipse from the references cited in || . It is important to remark here that although partial 
solar eclipses are not so useful to astronomers they are nevertheless relevant for neutrino 
physics as long as C is above 0.39. 

A study of the global maps of the five total/annular eclipses listed above shows that only 
three of them can occur as a single or a double eclipse at any of the three existing detector 
sites Kamioka, Sudbury and Gran Sasso with coordinates (36.4°iV, 140°i?), (46.5°iV, 81°H^) 
and (42.5°iV, 13.5°E) respectively. Approximate estimates are presented below for the du- 
ration and optical coverage of the relevant detectors at these sites: 

1997 March 9 - This eclipse will be a single eclipse for the Super-Kamioka, with an 
approximate duration of two and a half hours and a maximum optical coverage of 
just over 60%. In addition, there will be a double eclipse at the Gran-Sasso, the site 
of Borexino, at almost the same time as at Kamioka and for the same duration, with 
a 70% optical coverage. 

1998 Feb 26 - Though Sudbury will have an optical eclipse the coverage is less than 39% 
and so no neutrino-eclipse occurs. However, a double eclipse occurs at Gran-Sasso 
with a maximum coverage of 70% to 80%. This corresponds to the single eclipse at 
(31°A, 45°W). 

1999 Aug 11 - This eclipse will provide one of the best opportunities, as a single eclipse 
with 90% optical coverage and almost 3 hour duration at Gran-Sasso. In addition 
Super-Kamioka site will also get a double eclipse with 90% to 100% optical coverage, 
corresponding to the single eclipse at approximately (36° N, 40° E) 

III. THEORY 



5 



A. Regeneration in the moon 

We now describe a straightforward way of obtaining the neutrino regeneration effect in 
the moon by using a model of moon of constant density (3.33gms/cc). Let a neutrino of 
flavor a be produced at time t = to in the core of the sun. Its state vector is 

|^(to)) = |^)=E f/ f>f)- (6) 

i 

where |z/f) are the matter dependent mass eigenstates with mass eigenvalues /if and 
are the matrix elements of the matter dependent mixing matrix in the core of the sun. 
We use Greek index a to denote the three flavors e, /x,r and Latin index i to denote the 
mass eigenstates i = 1,2,3. The neutrino propagates in the sun adiabatically upto t R (the 
resonance point), makes non-adiabatic transitions at t R , propagates adiabatically upto t\ 
(the edge of the sun) and propagates as a free particle upto t 2 when it enters the moon. So 
the state vector at t 2 is 

= E Wi)exp (-iSj(t 2 - h)) exp (-i jf 1 e$(t)dt) M^exp (-i j* ef(t)dt) L/£. (7) 

where ef(t)(= E + (pf (t)) 2 /2E) are the matter dependent energy eigenvalues in the sun, Si 
and \vi) are the energy eigenvalues and the corresponding eigenstates in vacuum and Mjj 
is the probability amplitude for the non-adiabatic transition % — > j. We multiply the right 
hand side of eq.(^) by J2k l z/ fc f )( z/ fc^l ( = l)where (i =1,2,3) is the complete set of matter 
dependent mass eigenstates inside the moon. The neutrino propagates upto the the other 
end of the moon at t 3 , and the state vector at £3 is 
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We have introduced the probability amplitude M^j for non-adiabatic transitions j — > k due 
to the abrupt change in density when the neutrino enters the moon . It is given by 



= K M h-> = EK M K>KN = E u %v* ( 9 ) 

7 7 



where U^j is the mixing matrix in vacuum . We multiply the right hand side of eq. (|Sp by 
J2i KX^zl where \ui) (1 =1,2,3) is the complete set of vacuum mass eigenstates. The neutrino 
leaves the other end of the moon at t — t^ and propagates upto the surface of the earth 
,which it reaches at t± .So the state vector at t± is 

l*«(*4)> = E \"i) ex P (- i£ k(t3 ~ hj) (vi\vj?)Mj%exp (-isj(h -t x )-i t e^t)dt) 

h,j,i,l ^ JtR ' 



x M^exp ( -i jf R ef(t)dt ) U^exp ( - > :, ( t , - t,, } ) 



E HM^M^Mpiexpi-i^u) (10) 

k,j,i,l 



where 



= ef{h - t 2 ) + et{U - h) + Ej[t 2 - h) + / 1 (t)dt + [* ef{t)dt (11) 

Jtjj J to 



We have used the fact that the the probability amplitude for non-adiabatic transitions k — > / 
due to the abrupt change in density when the neutrino leaves the moon is 

faW> = «Uf (12) 
The probability of detecting a neutrino of flavor (3 at £4 is 

\(v0*{U))\ 2 = ^U^U^M^M^M^M^M^MU^exp(^ m - <% w ) (13) 

where the summation is over the set of indices k, I, k', I' Averaging over tR leads to 
Sa'Sjj' and this results in the desired incoherent mixture of mass eigenstates of neutrinos 
reaching the surface of the moon at ti- Calling this averaged probability as ( the 
probability for a neutrino produced in the sun as v a to be detected as up in the earth after 
passing through the moon), we can write the result as 
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P^ = J2 pS ( a ^j) pM U^(3) (14) 

j 

where 

^(a^J-) = EWlWl a (15) 
P M ti~+P) = E U^M^M^M^M^exp^-iie^-e^dM-^-e^r) (16) 

where we have replaced (i 3 — t 2 ) by g?m, the distance travelled by the neutrino inside the 
moon, and (£4 — £3) by r the distance travelled by the neutrino from the moon to the earth. 
If there is no moon , we put du = 0, so that P M (j — > (3) becomes \Upj\ 2 and so eqn.flUp 
reduces to the usual [PJT2] averaged probability for v a produced in the sun to be detected 
as vp in the earth : 

jS = EI^I 2 WW- (1 7 ) 

id 



B. Regeneration during double eclipse. 

We use a model of earth of constant density (5.52gms/cc). We start with ^ a (ti) given 
by eq. (pTOf) and multiply the right hand side by J2 P Wp)(Vp\ (— 1) where \v^) (i = 1, 2, 3) is 
a complete set of mass eigenstates inside the earth. The neutrino enters the earth at time t 
= i 4 and is detected at time t = t 5 inside the earth. The state vector at time t= t 5 is 

I*aft0> = E \y p )M^M^Mp s ai e W {-i^ klp ) (18) 

h,j,i,l,p 

where we have introduced the probability amplitude for non adiabatic transitions I — > p 
due to the abrupt change in density when the neutrino enters the earth . It is given by 

MS H<h> = EOC (1 9 ) 

and 

®ijki P = e*{h - U) + 4 f (£ 3 - £2) + £/(£ 4 - is) + ej(t 2 - h) + ef(t)dt + f* ef{t)dt (20) 
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The probability of detecting a neutrino of flavor (3 at £5 is 

x exp {-i{$ mp - $i>fk'i' P ') (21) 

where the summation is over the set of indices i,j,k,l,p,i',j',k',l'p' Again averaging over 
tfi and calling this averaged probability as P^p E ( the probability for a neutrino produced 
in the sun as v a to be detected as vp in the earth after passing through the moon and the 
earth), we can write the result as 

P$ E = Y.P 3 (<*-+3)P* tE (j-+0) (22) 
j 

where 

P ME (J - /?) = E Up*Up p , M E t M E *, My* M^l * My V 
i,k,p,i',k', P ' 

xexp (-i{e E - e E )d E - i(e^ - e#)d M - i(ej - £z') r ) (23) 

where we have replaced (t 5 — t 4 ) by d E , the distance travelled by the neutrino inside the 
earth, (£3 — £2) by d,M, the distance travelled by the neutrino inside the moon and (£4 — £3) 
by r the distance travelled by the neutrino from the moon to the earth. 

For the sake of completeness, we state that if we put du = in eqn(p3|) we get the 
regeneration in earth alone: 



P*p = i:P S (a^3)P E ti^P) (24) 
j 

where 

P E U - = E UftUfcM%M%exp {-i{e E - e E )d E ) (25) 
k,k' 

Eqns([?|) and (|25| ) have been used to study the day-night effect 0. 



It is important to note that the factorization of probabilities seen in eqs(|HD,(22"D and 
is valid only for mass eigenstates in the intermediate state. An equivalent statement of 
this result is that the density matrix is diagonal only in the mass-eigenstate representation 
and not in the flavour representation. 



C. Three flavor mixing parameters 



We parameterize the mixing matrix U in vacuum as U = U 23 (ip)U 13 (<p)U 12 (u) where 
U^(8ij) is the two flavor mixing matrix between the i th and the j th mass eigenstates with 
the mixing angle 8{j, neglecting CP violation. In the solar neutrino problem if) drops out 
I~0| , [TT|] The mass differences in vacuum are defined as 5 2 i = //| — /xf and S 3 i = fi\ — fif. It 



has been shown |[T^JT^| that the simultaneous solution of both the solar and the atmospheric 
neutrino problems requires 

631 > 621 (26) 

and under this condition $31 also drops out. The rediagonalization of the mass matrix in the 
presence of matter (in solar core, moon or earth) under condition (p6| ) leads to the following 
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S21 sin 2u 
S21 cos Ioj — A cos 2 
sin0 m = sin0 (28) 



tan2w m = 2 - — (27) 



$21 = ^21 cos 2(cu — uj m ) — A cos 2 (j) cos 2to m (29) 

where A is the Wolfenstein term A = 2a/2 Gp N e E (N e is the number density of electrons 
and E is the neutrino energy) . We note that 8 31 > A S ,A M ,A E . In eqs.([Z7|)-(|29|) the 
subscript "m" stands for matter. Under the condition 631 ^> A w 5 2 i we need the non 
adiabatic transition probability |M^| 2 for i,j =1,2 only and this is taken to be the modified 
Landau-Zener jump probability for an exponentially varying solar density flU]] . 

IV. CALCULATIONS AND RESULTS 

The neutrino detection rates for a Super-Kamoika type of detector is given by 

R = J (j){E) a(E)P ee dE + ~(| (<j)(E)a(E)(l - P ee )dE) (30) 
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where the second term is the neutral current contribution and (j)(E) is the solar neutrino 
flux as a function of the neutrino energy E and cr(E) is the cross section from neutrino 
electron scattering and we integrate from 5MeV onwards. The cross section is taken from 
fl~4| and the flux from JOJ . The rates for a single eclipse, double eclipse and without eclipse 
(at day-time) Rm, Rme and Ro are calculated using Pjf, P^ E and P° e from eqns.(|T4"|),([2"iZD 
and ([17]) respectively. We define the enhancement factors F and G for a single and double 
eclipse respectively: 

F = Rm ~ R ° (31) 
Ro 

G = Rme ~ R ° . (32) 
Ro 

It is easy to see that F and G have to be less than 5 and this theoretical maximum value 
occurs when P° e = and Pjjf and P^ E are put 1. If one imposes the constraint that the 
observed |I(J neutrino rate is 0.51 ± 0.07 of the prediction of the standard solar model, the 



maximum possible enhancement is reduced to about 1.40 (at 90% C.L.). 

We calculate the enhancement factors F and G for various values of the neutrino param- 
eters, u, 621, and 0. We show the results as contour plots in the 821-uj plane for different 
values of 0. Figs|3] and |] show the F-contours for = 0° and = 30° respectively. For each 
we show the contours for different distances of travel of the neutrino through the moon. 
Fig. [| shows the G-contours for = for the maximum distance of travel of the neutrino 
inside the moon and the earth. The main features of the results are as follows: 

• As the distance travelled by the neutrino inside the moon increases one can see an 
appreciable increase in the enhancement factor F. It increases from less than 10% to 
about almost 100% when the neutrino travels the whole diameter of the moon in the 
case of two flavor mixing i.e = 0. 

• Large (> 40%) values of F occur for 10 between 20° and 30° and b?\ ~ 10~ 6 el/ 2 This 
is true even for nonzero 0. 
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• The effect of a non zero "13" mixing angle 0, is to dilute the enhancement factor F for 
all values of distance travelled through the moon. (In fact for ss 45°, F is practically 
zero and so we have not plotted this case.) This is because a non zero means u e <-> u T 
oscillations, and matter cannot reconvert v T back to v e , because the "13" mixing angle 
is not affected by matter. 

• If large enhancement F is seen for values of x < 0.6, it immediately signals a very 
small value of 0. On the other hand, if no enhancement is seen for small x but there 
is enhancement only for x > 0.8 it signals an appreciable value of 0. 

• For a double eclipse there are considerable enhancements even for small values of uj. 
There is enhancement throughout the range of uj from small angles till about 40°. In 
fact the regions of largest enhancement (> 100%) are for uj between 5° to 20°. 

• The region of maximum enhancement factor G is centered around a value of 5 2 i which 
is a little above lO -6 ^ 2 , this being the value for maximum in F. This can be traced 
to the fact : A E > A M . However sizeable enhancement occurs over a wide range of 

$21- 

• If enhancement is not seen, then certain regions of the neutrino parameter space can 
be excluded. If no enhancement is seen for single eclipse, a panel of uj between 5 — 25° 
and 5 2 i ~ 2 x 10 -7 — 2 x lO^eV 2 for = can be ruled out. If it is not seen for a 
double eclipse, a larger region can be ruled out. 

V. DISCUSSION 

We have studied the effect on the solar neutrinos of their passage through the moon as 
well as the moon together with the earth. Although the numerical results presented in the 
paper cover only a representative sample of the set of various parameters, our analytical 
expressions can be used for more extensive calculations depending on the requirement. Also 
one can go beyond the hierarchy : 5 3i >> 5 2 i- 
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We now offer a few concluding remarks: 

• Together with the day-night effect, the eclipse effects provide us with the tools for 
studying solar neutrinos, in a way independent of the uncertainties of the solar models. 

• If the neutrino mass differences are really very small (S21 < 10~ 5 eV 2 ) there is no way 
of pinning down the neutrino parameters other than using the astronomical objects 
such as the moon or the earth for the " long-base-line experiments" . 

• It is important to stress that even the demonstrated absence of any eclipse effect would 
provide us with definitive information on neutrino physics. 

• Accumulation of data over many eclipses may be needed for good statistics. In any 
future planning of detector sites, this may be kept in mind. 

• It appears that Nature has chosen the neutrino parameters in such a way that the sun 
affects the propagation of solar neutrinos. It may be hoped that Nature has similarly 
chosen "lucky" parameters so that the moon and the earth too can affect the neutrinos! 

• Finally, we stress the novelty of the whole phenomenon, and urge the experimentalists 
to look for and study the eclipse effects in an unbiased manner. They may even discover 
some surprises, not predicted by our calculations! 

Acknowledgments: We thank KVL Sarma for bringing ref to our attention, N.D. 
Hari Dass for raising the possibility of detecting the neutrinos on the opposite side of the 
earth during the eclipse, and M.C Sinha for much help with the astronomical aspects of the 
problem and Sandip Pakvasa for useful discussions. 
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FIGURES 

FIG. 1. The fractional distance travelled by the neutrino inside the moon x(= ) is plotted 
against the optical coverage C of the solar eclipse. 

FIG. 2. Geometry relating the double eclipse point (A, a) to the single eclipse point (a, (3). (a) 
Section of the earth passing through (a, (3) and perpendicular to the ecliptic, (b) Section passing 
through (a, (3) and parallel to the equator. 

FIG. 3. Contour plots of the enhancement factor for single eclipse F(= ^ M ) in the 

Ro 

u — 821 plane for = 0° and for four values of x (x = 0.4, 0.6, 0.8 and 1.0). The enhancement 
factor (regarded as a percentage) increases by 10% for every adjacent ring , as we move inwards 
towards the centre of the plot. 

FIG. 4. Contour plots of the enhancement factor for single eclipse F{= — ^— ) in the 

Ro 

w ~ ^21 plane for <fi = 30° and x = 0.6, 0.8 and 1.0. The enhancement factor (regarded as a 
percentage) increases by 10% for every adjacent ring , as we move inwards towards the centre of 
the plot. 

FIG. 5. Contour plot of the enhancement factor for double eclipse G{= — ——^ — ) for <\> = 

Ro 

and x= 1.0. The distance travelled by the neutrino inside the earth is also taken to be the full 
earth diameter. The enhancement factor increases by 20% as we move inwards. 
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